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Abstract. In this talk we present results for the propagators and the phase diagram of 
ON . QCD, obtained from a truncated set of Dyson- Schwinger equations at finite temperature T 

and quark chemical potential /j,. We include back-coupling effects of the dressed quark and 
gluon propagators, which also allows us to study the influence of strange quarks. For the phase 
' (— i ', diagram we find a critical end-point at fi/T ~ 1.9 and coinciding chiral and deconfinement phase 

Q • transitions. 

CD 

1. Introduction 

The most reliable source of information on quantum chromodynamics (QCD) at finite 
temperature has been lattice QCD which finds a first order phase transition in the quenched 
t^- ■ case and a crossover when dynamical quarks are included [HE]- hi the latter case the order 

parameter for confinement, the Polyakov loop, shows a crossover transition nearby the chiral 
phase transition. However, at finite density lattice studies are hindered by the fermion sign 
problem. This led to a large industry of studies in effective models such as the Polyakov-loop 
; extended Nambu-Jona-Lasinio [21 H] and quark-meson models [SJ El E] • From these studies the 

general expectation that the crossover turns into a first order phase transition at a critical end- 
point (CEP) has emerged, with the position and even the existence of the CEP under debate. 
When employing effective models one has to trade the QCD degrees of freedom for effective ones. 
In contrast to that we use Dyson-Schwinger equations (DSEs) to study QCD and its degrees 
of freedom, quarks and gluons, directly. We take advantage of this possibility by introducing a 
truncation where the back-reaction of quarks onto gluons is explicitly taken into account. Here 
we summarize our results published in Ref. [14]. 

2. Order parameters from Dyson-Schwinger equations 

The main objects of interest in this work are the propagators for the fully dressed quark and 
gluon, which are given, respectively, as 

S(p) = [i(uj n + ifi)j 4 C(p) + ip^A(p) + B( P )]- 1 , 

(1) 

DAP) = P»^ + P»^^ (2) 

where p = {p, u p ). The quark is dressed by the scalar functions A, B and C which depend on 
and Up only. The in-medium gluon has two dressing functions Zl and Zt for the components 
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Figure 1. The coupled Dyson-Schwinger equations for the quark and gluon propagators. 

longitudinal and transversal to the medium. P^ u and Pj^ u are the corresponding projectors. In 
the quark propagator we neglected a fourth dressing functions which comes with jryj^ and can 
be shown by explicit calculations to have virtually no influence on the results. 

2.1. Quark condensate 

In the quark propagator, Eq. ([1]), the presence of chiral symmetry breaking is directly reflected 
in a non- vanishing scalar dressing funtion B. A derived order parameter for chiral symmetry 
breaking is the condensate 



which is calculated from the propagator once the quark DSE is solved. 
2.2. Dressed Polyakov loop 

It is a much more non-trivial task to find an order parameter for confinement. The desirable 
object would be the Polyakov loop which can be connected to the free energy of a static quark, 
and is thus a description for confinement. Since we can not directly access the Polyakov loop 
from solutions of DSEs we will use the so-called dual condensates. They are defined as 



where (ipip)^ is the quark condensate evaluated at generalized £/(l)-valued boundary conditions 
ip(x,l/T) = e i( ^ijj(x,0) with cp € [0,2vr[. In [SJ El ED] it has been shown that £ ±1 contains 
the ordinary Polyakov loop as well as contributions from spatial detours. It has therefore been 
called the dressed Polyakov loop. Most notably it is sensible to center symmetry, just like the 
conventional Polyakov loop. It is therefore an order parameter for confinement that is accessible 
by continuum functional methods, since only the dressed quark propagator is necessary as an 
input [TT] . 

With Eqs. (|3l4p we have identified order parameters that rely only on the quark propagator. 
We will now define the truncation scheme that we will use to get our hands on not only the 
quark but also the gluon propagator. 
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3. Truncation scheme 

The upper diagram in Fig. (pQ) shows the quark DSE. This equation defines the fully dressed 
quark propagator. On the right hand side we have the self-energy which depends on the fully 
dressed gluon and quark- gluon vertex. We need a prescription for these objects in order to 
get a closed set of equations, this will be given in the following. The challenge of in-medium 
Dyson-Schwinger equations is to include the relevant medium effects in the applied truncation 
scheme. We will here use the same construction for the quark-gluon vertex as in \1'2\ [T3 1 114j: 



here q is the gluon momentum while p and k are the quark momenta. This ansatz contains the 
first part of the Ball-Chiu construction and is thus dependent on the medium. On the other 
hand, the dressing function T is constant in T and /i. This works well up to T c but leads to an 
over-estimation of chiral symmetry breaking in the quark-gluon plasma as will be shown below. 

The parameters are fixed to reproduce the pion decay constant in the vacuum and turn out 
to be di = 7.5 GeV 2 , d 2 = 0.5 GeV 2 and A = 1.4 GeV. In the UV we have <5 = -9 44N ff 8Nf and 

Q UN c -2N f 

For the gluon propagator our strategy is to start with quenched QCD. Here, lattice 
calculations are up to now the most reliable source for the Landau gauge and fully temperature 
dependent propagator [15]. We have used this input to calculate the first (second) order phase 
transition in SU(S) (SU(2)) gauge theory in [IB] . To be able to correctly describe the crossover 
in unquenched QCD and to go to finite chemical potential, we add the quark loop from the 
gluon DSE to the quenched propagator. In Fig. ([T]) the resulting set of equations that we have 
solved is shown diagrammatically, where the first term on the right hand side of the gluon DSE 
denotes the quenched propagator. A strange-quark loop is included for the Nf = 2+1 case. This 
procedure for unquenching the gluon propagator neglects all quark loops inside the Yang-Mills 
self-energies. In the vacuum this leads to an error on the 5 percent level. 

The quark loop at finite temperature can be split into longitudinal (L) and transversal (T) 
parts: 

rv(p) = pT v ii T (p) + p^n L ( P ), (7) 

where one notes that lim p ^.Q p 2 H T (p) = but 

lim/n L (p) = ^ (8) 
P^O 2 

i.e. the longitudinal part produces a thermal (or Debye) screening mass. In the color 
superconducting phase the transversal part would generate a Meissner mass in a similar way. 

4. Results 

We present results from a calculation of the quark condensates and the thermal mass for 
N f = 2 + 1 QCD in the left side of Fig. Q. 

What we find from the light quark condensate is a crossover for small chemical potential, 
which becomes stronger with rising \i and eventually turns into a first order phase transition. 
The strange quark condensate reflects this change due to the coupling to the light quark via the 
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Figure 2. The left figure shows the light and strange quark condensates as the lower and upper 
surfaces, respectively. The right figure shows the thermal mass contribution from the quark 
loop, normalized by T 2 + §/x 2 , the asymptotic value. Both figures are adopted from Ref. [14j . 

gluon propagator. However, it continues to decrease at larger temperatures due to the larger 
strange quark mass. 

The contribution to the thermal gluon mass due to the quark loop, Eq. ([8]), is shown to 
the right in Fig. , normalized by its asymptotic HTL/HDL behaviour. It is small in the 
hadronic phase, where it is suppressed by the inverse quark mass. Around the phase transition 
it rises strongly and approaches the asymptotic value above T c . For large /i, at the first order 
phase transition, we observe a jump, which is inherited from the jump in the quark propagator. 
The quark loop reduces the gluon strength, and therefore also the quark self-energy. A smaller 
self energy in turn leads to a larger quark loop, and therefore the back-coupling of quark and 
gluon has an accelerating effect on the phase transition. This acceleration will be important to 
reproduce lattice data, which we do in the following. 

From the light and strange quark condensates we can obtain the quantity 

A ltS = " — {W)„ (9) 

m s 

which is finite since the divergent term mA 2 in ('(pip) cancels. It can therefore be used in 
comparisons with different approaches like lattice QCD. This comparison is shown in the left 
part of Fig. Q. It shows a rather good agreement up to T c and an overestimation of the 
remaining chiral symmetry breaking in our approach above T c . This we attribute to the constant 
vertex strength, which should be temperature dependent in a refined model. Nonetheless we 
find a critical temperature of T c ~ 156 MeV in good agreement with the lattice pQ, and a similar 
steepness of the crossover. 

Finally we show the main result of Ref. |14j . the phase diagram of Nf = 2 and Nf = 2 + 1 
QCD in the right part of Fig. Q. We find a critical end-point at 0, T) « (190 MeV, 100 MeV) 
for Nf = 2 + 1, which is certainly outside the area that can be accessed by extrapolations from 
lattice data. The (pseudo-)critical temperature for deconfinement lies slightly above that of 
chiral restoration in the crossover regime, and coincides close to and at the CEP. Fig. ([3|) also 
shows in lighter colors the phase diagram for the two flavor case, where T c ~ 200 MeV. The 
impact of the strange quark is evidently a reduction of T c by about 50 MeV while the CEP 
moves only slightly to larger [i. 




Figure 3. In the left figure we compare the regulated and normalized quark condensate to 
corresponding results from lattice QCD [I]. In the right figure we show the phase diagram for 
Nf = 2 + 1 QCD. In lighter colors we also show the results for Nf = 2 QCD (upper lines). Both 
figures are adopted from Ref. |14j . 

5. Summary 

We discussed a truncation scheme to the Dyson-Schwinger equations of QCD that explicitly 
couples the quark and gluon propagators. This makes the gluon sensitive to the chiral phase 
transition via the thermal mass contribution from the quark loop. In this approximation we 
showed that the condensate is in reasonable agreement with lattice simulations. For the phase 
diagram we find a CEP at fj,/T ~ 1.9, and nearby chiral and deconfinement phase transitions. 
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